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We demonstrate a "bottom up" approach to the com- 
putational design of a multifunctional chemical sensor. 
General techniques are employed for describing the ad- 
sorption coverage and resistance properties of the sensor 
based on density functional theory (DFT) and non-equi- 



librium Green's function methodologies (NEGF), re- 
spectively. Specifically, we show how Ni and Cu doped 
metallic (6,6) single-walled carbon nanotubes (SWNTs) 
may work as effective multifunctional sensors for both 
CO and NH3. 



1 Introduction Detecting specific chemical species 
at small concentrations is of fundamental importance for 
many industrial/scientific processes, medical applications, 
and environmental monitoring [1]. Nano structured ma- 
terials are ideally suited for sensor applications because 
of their large surface to volume ratio, making them sen- 
sitive to the adsorption of individual molecules. Specifi- 
cally, single-walled carbon nanotubes (SWNTs) [2] work 
remarkably well as detectors for small gas molecules, 
as demonstrated for both individual SWNTs [3-8] and 
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Figure 1 Schematic of a chemical sensor consisting of ac- 
tive sites (metal dopants in a (6,6) carbon nanotube), a tar- 
get molecule (CO), a background (atmospheric air), and a 
sensing property (resistance). 



SWNT networks [9, 10]. Previous studies have shown 
that SWNTs are highly sensitive to most molecules upon 
functionalization [11-19]. However, the difficulty is deter- 
mining which specific molecules are present. In this study 
we show how changing the functionalization provides "an- 
other handle" for differentiating the SWNTs response of 
different gases/molecules. 

Recent experimental advances now make the con- 
trolled doping of chirality selected SWNTs with metal 
atoms a possibility. Specifically, these include (1) photo- 
moluminesence, Raman and XAS techniques for measur- 
ing the fraction of various SWNT chiralities in an enriched 
sample [20-22]; (2) the separation of SWNT samples by 
chirality using DNA wrapping [23-26], chromatographic 
separation [26, 27] and Density Gradient Ultracentrifu- 
gation (DGU)[28]; (3) SWNT resonators for measuring 
individual atoms of a metal vapour which adsorb on a 
SWNT [29, 30]; and (4) aberration corrected low energy 
(<50 keV) transmission electron microscopy (TEM) [31]. 
The latter provides control over the formation of defects 
in situ by adjusting the energy of the electron beam above 
and below the threshold energy for defect formation at 
a specific location in the SWNT. These methods provide 
such a high level of control that it is now possible for ex- 
perimentalists to take a specific SWNT chirality and dope 
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Figure 2 Binding energies in eV for N2, O2, H2O, CO, 
NH3, and H2S on a Ni atom occupying a divacancy site 
(grey) and on a Cu atom occupying a monovacancy site 
(orange) in a (6,6) carbon nanotube. Schematics of the two 
active sites are shown above. 
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Figure 3 Fractional coverage in thermal equilibrium of 
Ni in a divacancy (top) and Cu in a monovacancy (bot- 
tom) versus CO concentration (left) and NH3 concentra- 
tion (right) in a background of air at room temperature and 
1 bar of pressure. 



the structure with individual metal atoms at a specified 
location. 

At the same time, theorists are now able to embrace 
a "bottom up" approach to the design of nanosensors, 
harnessing the thermodynamics of self-assembly to find 
useful sensing systems in silico. With recent advances in 
both computational power and methodologies, theorists 
can now efficiently and accurately screen hundreds of 
candidate sensor designs using a combination of density 
functional theory (DFT) for energetics of adsorption and 
stability, and non-equilibrium Green's function (NEGF) 
methodologies for the electrical response [32]. 

In general, any nanosensing system consists of the fol- 
lowing four main components: (1) a "target molecule" to 
be detected, (2) an "active site" where the target molecule 
may adsorb on the sensor, (3) a property of the sensor 
which changes upon adsorption of the target molecule, and 
(4) a "background" of adsorbing molecules which make 
up the background signal, as depicted in Fig. 1. For the 
sensor to be effective, the active site must be designed so 
that adsorption of the target molecule in the presence of the 
background is sufficient to change the sensing property. 

In this study we shall consider Ni and Cu doped metal- 
lic (6,6) SWNTs as our active sites, specifically when Ni 
replaces two C atoms in a divacancy structure, and Cu re- 
places a single C atoms in the monovacancy structure, as 
shown schematically in Fig. 2. Of the two possible ori- 
entations for a divacancy in an armchair SWNT, we here 
consider the one in which Ni adsorption has been shown 
to be more favourable. These two active sites have pre- 
viously been shown to be both stable, and somewhat re- 
sistant to oxidation under standard atmospheric pressure 



Table 1 Equilibrium atmospheric concentrations at 4% ab- 
solute humidity. 



Background 


Concentration 


Target 


Concentration 


N2 


74.96% 


CO 


96.00 ppb 


O2 


20.11% 


NH3 


16.32 ppb 


H2O 


4.00% 


H2S 


0.96 ppb 



and room temperature [32]. As such, we shall consider the 
background to have a standard atmospheric composition 
with an absolute humidity of 4%, shown in Table 1. 

2 Methodology All DFT calculations have employed 
the PBE exchange correlation (xc)-functional [33]. Ad- 
sorption energies and structural minimization calcula- 
tions have been performed using the real-space projector- 
augmented- wave method DFT code GPAW [34], with a 
grid spacing of 0.2 A, taking spin polarization into ac- 
count. NEGF transport calculations [35] have employed 
the electronic Hamiltonian obtained from the DFT code 
SIESTA [36], with a double zeta polarized (DZP) locaHzed 
basis set. 

Vacancies in the (6,6) metallic carbon nanotube have 
been modelled using six repeated minimal unit cells, with 
a supercell oflSAx ISAx 14.622 A. Such a large super- 
cell was necessary to both minimize interactions between 
defects and adsorbates and to ensure the Hamiltonian of the 
nanotube layers adjacent to the cell boundaries was con- 
verged to within 0.1 eV of that for the pristine nanotube. 
For such a large supercell, a F point calculation was found 
to be sufficient to describe the nanotube 's periodicity. 
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(a) Cu in Monovacancy (b) Ni in Divacancy 
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Figure 4 Density of states (top) and transmission function 
(bottom) versus energy e relative to the Fermi level in 
eV for clean (solid line; black), O2 occupied (dashed line; 
red), NH3 occupied (dash-dotted line; orange), and CO oc- 
cupied (dotted line; green) (a) Cu in monovacancy and (b) 
Ni in divacancy II of a (6,6) carbon nanotube. 
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Figure 5 Change in resistance AR in Q (top) and mfl 
(bottom) per active site versus CO (left) and NH3 (right) 
concentration in ppm for Ni in a divacancy II (dashed line) 
and Cu in a monovacancy (solid line) of a (6,6) carbon nan- 
otube. The reference concentrations of CO and NH3 are 0. 1 
and 0.01 ppm, respectively, in a background of air at room 
temperature and 1 bar of pressure. 



Binding energies £^bind are defined as the difference 
in total energies between a gas phase target/background 
molecule with a clean active site and a target/background 
molecule adsorbed on an active site. As seen in Fig. 2, CO 
has a strong binding energy to both active sites, while NH3 
binds only to Cu in a monovacancy site. We shall use this 
in our design of a multifunctional sensor for both CO and 
NH3. 

3 Results & Discussion For a sensor containing 
hundreds, if not thousands of sites, the fractional coverage 
of sites is reasonably described by the fractional coverage 
in thermodynamic equilibrium O at standard temperature 
and pressure in terms of the target molecule concentra- 
tion. Specifically for a sufficiently large ensemble of active 
sites, the fractional coverage of sites by a given adsorbate 
X in a background set of molecules B at equilibrium is 
given by 

K[X]C[X] 



0[X] 



(1) 



EYeBK[Y]C[Yy 

where K = fe+//c_ is the ratio of forward and backward 
rate constants for the adsorption reaction. This may be ex- 
pressed as 

^bind [X] — TSgas [X] 



K[X] = exp 



kr,T 



(2) 



where 6'gas[-^] is the gas phase entropy of species X [37], 
/cb is Boltzmann's constant and T is the temperature. 

In Fig. 3 we see the coverages of both Ni and Cu active 
sites as a function of CO and NH3 concentration. Clearly, 



both Cu and Ni active sites will be sensitive to the adsorp- 
tion of CO. On the other hand, the coverage of the Cu ac- 
tive sites is sensitive to NH3 concentration, while Ni active 
sites are not. This suggests that by combining the response 
of Cu and Ni active sites, we may obtain a multifunctional 
sensor for both CO and NH3 . 

However, although molecules may adsorb on the active 
sites, the sensing property of the SWNT resistance must 
change significantly for the sensor to be useful. In Fig. 4 
we show both the density of states (DOS) and transmission 
probability T(s) for an electron of a given energy e to travel 
past an active site as calculated within the NEGF method. 
Since the energy levels of a given molecule may be thought 
of as its electronic "finger print", so long as there is suf- 
ficient binding between the molecule and active site, this 
"finger print" is evidenced in the DOS of the system. Fur- 
thermore, the peaks of the DOS due to the molecule tend 
to form Fano anti-resonances in the transmission, as elec- 
trons are scattered off these states (c.f. Fig. 4). In this way, 
any adsorbed molecule tends to leave its "finger print" 
in the transmission probability through an active site. For 
this reason, conductance/resistance is generally an effec- 
tive sensing property. 

Finally, we may estimate the average change in resis- 
tance AR of mi active site as a function of target molecule 
concentration. As discussed in Ref. [32], this change in re- 
sistance is reasonably well described by 



AR^Y. Rs{x){e[x, c] - e[x, Co]), 



(3) 



X 
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where C is the concentration, Co is the concentra- 
tion at standard temperature and pressure, Rs{X) = 
Go{l/T{eF,X) - 1/(2), Go = 2e^/h is the quantum 
of conductance, and T{eF,X) is the transmission proba- 
biUty at the Fermi level through an active site with species 
X adsorbed. 

In Fig. 5 we show AR for a single active site (Ni or 
Cu) as a function of target molecule concentration (CO 
or NH3). Keeping in mind that a change in concentration 
from 1 ppm to 50 ppm amounts to a change from allowable 
to toxic concentrations for both CO and NH3 , we see that 
this sensor design should work effectively for both target 
molecules. 

4 Conclusions In conclusion, we have demonstrated 
in silico how a combined Cu and Ni doped metallic SWNT 
device may work effectively as a multifunctional sensor for 
both CO and NH3. The methodology employed is appli- 
cable to other sensor designs and different environments. 
Also, by varying the metal dopant, we obtain "another han- 
dle" for tuning the sensitivity of our sensor. 
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